Unitary expansion of the time evolution operator 
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We propose an expansion of the unitary evolution operator, associated to a given Schrodinger 
equation, in terms of a finite product of explicit unitary operators. In this manner, this unitary 
expansion can be truncated at the desired level of approximation, as shown in the given examples. 
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I. INTRODUCTION 

The central problem in any dynamical theory is to find 
the time evolution of a system that was prepared in a 
given initial state. In quantum mechanics there are only 
a few of these problems that are readily solved by simple 
analytical methods [1 . In general, we have to rely on 
approximations to obtain out of the Schrodinger equa- 
tion the time evolution operator Ug(t,to) in a suitable 
form. With the explicit knowledge of Us(t,t ), we may 
calculate the expectation value of any physical observ- 
able of the system at any time t once we know the state 
of system at the time to- Frequently, we are able to find 
the time evolution operator Uo(t, to) associated with Hq, 
a part of the total Hamiltonian H = Ho + Hi nt . In this 
case, it is usually convenient to make a transformation 
to the "interaction picture" such that 



U s (t,t ) = U (t,to)Uj{t,to) 



(1) 



holds with Uo(to,t ) = Ui(t n ,t n ) = I. Consequently, 
the time evolution operator in the interaction picture, 
Ui(t,to), satisfies the Schrodinger equation 



dUi(t,t Q ) 

at 



-i\H x {t)Ui{t,t Q ) 



(2) 



where we have considered h — 1 and have defined 
Afli(i) = Ul(t,to)H int (t)Uo(t,t ). The parameter A is 
some dimensionless real parameter chosen to give a mea- 
sure of the relative magnitude of the most important ma- 
trix elements of Hi nt and Hq in a given problem. 

The most popular perturbation approximation method 
to deal with the Schrodinger equation in Eq. ^ is the 
Dyson expansion [2]: 

Ui(t,t ) = T-i\ [ dtiHxih, t ) + 
Jo 

HA) 2 f dh f dt 2 H 1 (t 1 ,t Q )H 1 {t 2 ,t Q ) + (3) 
Jo Jo 

where it is very convenient to estimate the solution by 
truncating the series to a given order A™. Besides the nor- 
mal difficulties to calculate high-order terms, the Dyson 
truncation produces an approximated evolution operator 



that is not unitary. Other expansions of operator U\{t, to) 
have also been proposed in the literature, as the Magnus 
expansion [3] the Fer product [3] and, more recently, the 
Anicllo expansion [5]. New bounds for the convergence 
of the Magnus expansion and the Fer product have been 
recently studied in Ref. [6]. Other product expansions 
have also been considered in the literature [7J. 

In this paper, we present an alternative method to cal- 
culate the time evolution operator Uj(t,to) as a product 
of a finite number of unitary operators 

Uj(t, t ) = Ui(t, t )U 2 (t, t )--- U k {t, t )--- U N (t, t ) , 

where each operator Uk(t,to),k — 1,2, ...N — 1, can be 
written as the exponential of an anti-Hermitian opera- 
tor proportional to (A) fc , while J7jy(t, to) — ^ is at most 
of order (A)^. The number N of operators in the ex- 
pansion can be as large as we want. If we approximate 
the last operator Uj^(t,to) in the product by the unit 
operator, we obtain an expansion of the evolution opera- 
tor Ui(t,to) which is explicitly unitary to order (A)^ 1 . 
Besides this important advantage, this expansion is well 
suited to treat a variety of problems. In Section II, we 
derive the expressions for each term in the expansion; 
in Section III, we provide pedagogical examples; and in 
Section IV, we present our conclusions. 



II. THE METHOD 

We start with the simple case of N — 2. Equation Q 
can then be written as 



[//(Mo) = Ui(t,to)U 3 (t,to). 



(5) 



Following the same kind of procedure used in the trans- 
formation to the interaction representation, we write 



U I (t,t ) = e- iXW ^U 2 (t,t Q ) 



(6) 



where W\{t,to) is an hermitian operator to be chosen 
conveniently. From now on, we set t = and avoid 
writing it to simplify the notation. From Eqs. ([2]), ^ 
and ([6]), we have 



dU 2 {t) 
dt 



-i\H 2 {t)U 2 {t) , 



(7) 



2 



where 



ff 2 (*)=x; i M-(« i ^(*)) m {^(*) 

'-^ ml I 



1 dWi(t)- 



m + 1 e>t 



(8) 



Here, we have defined <2gL4{-} = [A,-] and used the fol- 
lowing relation 



e aA 5e- Qi = E ^{adA) m {B} 



B + a[A,B] + —[A,[A,B]] + ~. (9) 



Choosing 



we have 



Wi(t) = / Hx(t')d£ 
Jo 



(10) 



^(*) - E A^M^(t)r{mffiW} , (ii) 

m— 1 ^ '* 

which is of order A. In this case, [^(t) is the solution 
of Eq. ([7]) and should be written as an exponential of a 
non-Hermitian operator that is, in general, a series on 
the variable A, starting with A 2 . 

In the simple case where [Hx(t), H\(f)\ — —2if(t,t') 
is a c-number function, then 



H 2 (t) = X / dt'f{t',t) 



(12) 



is also a c-number function. Consequently, Eq. ^ can 
be easily integrated to give U%(t) = e~ lX where 



0(i)= / / dt"f(t",t') 
Jo Jo 

and the time evolution operator is 



Uj(t) = e" iA -^ d *i^i(ti) e -^ 2 <? i ( t ). 



(13) 



(14) 



This result is well known and could also be easily ob- 
tained using the Magnus expansion [8]. It can be used, 



for example, to easily obtain the time evolution operator 
for the quantum state generated by an external time- 
dependent force acting on a mechanical oscillator, as we 
show in the next section. 

The procedure described above can be generalized for 
any value of N greater than 2 by setting 

U n (t) = e- lX ^Mu n+1 (t) , n = 2,3..JV - 1 , (15) 

so that expansion of the operator Ui(t) may also read 

Uj(t) = e -i*W e -*W... e -i*- 1 (^ JV (t). (16) 
The operators U n (t), for n = 2, 3.. satisfy a Schrodinger- 
like equation 



dU n (t) 

dt 

where H n (t) is given by 
(iX) 



-iXH n (t)U n (t) 



(17) 



Hn(t) = J2 



. {adW n -i(t)) r 
ml 



:{fl»-i(t) 



1 dWn-i(t) 



1 



dt- 



} ■ (18) 



By choosing operators Wj(t)'s for j = 1, ...N — 1, we 
obtain operators Hj(t) for j = 2, ...N, and the expansion 
given by Eq. (16). 



We now show that it is possible to choose opera- 
tors W n (t) being proportional to A™ -1 , and such that 
the operators H n (t) are power series in the variable A 
starting with the power A" -1 . Then, by substituting 
W n {t),n = 1,2. ..(TV - 1) in Eq. Q and noticing that 
/ — U^(t) would be at least of (9(A ), we will obtain the 
desired expansion announced in the Introduction. We 
make the proof by construction. Writing explicitly the 
dependence of the operators W n (t) and H k (t) ona = iX, 
we have 



W k {t) = op- y W k {t) , 

OO 

H k (t) = E H kJ (t)c 

j=k-l 



(19) 



By substituting Eq. ( 19 1 in Eq. ( 18 ) we get, for j > n > 1, 



OO Tb OO -4 IV -| 

m=lfe=li=/s-l ' fe=lm=0 ^ 

and for < j < n, 



5(fem+fc— 1— j)\ , 



(20) 



dt 



OO 71 OO - ib uu - _JTT7" 

EE E ^w(*)r{i^^ 

m-l fc-1 k — 1 m— 1 

(21) 



3 



where 5(m — n) = S m ^ n i s the Kronecker delta. Notice 
that j + 1>A; in Eq. ^ so that dW ^ t l{t) is given re- 
cursively in terms of the operators Wk(t) and Hk(t), for 
k < j. For example, if we set n = 2 and j = 1 in the 
above equations, we easily get 



dW 2 {t) 
dt 



where we used Eq. ( 21 ) and the fact that Hifi(t) 



Using the initial condition H^(0) = 0, we have 



W 2 {t) 



1 



2 jo 

which also can be written as 

W 2 (t) = \ f dh / 1 dt 2 [^i(*2),ffi(*i)] 



(22) 

(23) 
(24) 



To obtain an approximate expressio n for Ui(t) valid to 
order 0(A 2 ), we first set N = 3 in Eq (fl6|: 



tTr(f) = e - lXW ^e-^ 2w ^Uz{t) . (25) 
f/ 3 (t) — I is of order A 3 , since it satisfies the Schrddinger 



equation, Eq. (17), with H 3 (t) of the order 0(A 2 ). If we 



approximate U 3 (t) by the identity we get an approxima- 
tion which is unitary and valid to order C(A 2 ). Using the 

expressions for W\(t) = W\(t) and for W2(t) given in 
Eq. flOl) and Eq. (l23|), we have 



U^t) « exp{(-iA) / dtiHi{t x )} x 
Jo 

exp{— / dii / dt 2 [H 1 (t 2 ),H 1 (t 1 )]}. 
z Jo Jo 



(26) 



The procedure described above can be generalized for 
obtaining approximations involving a product of N op- 
erators, by calculating Wk(t),k — l,...N, through Eqs. 
(20) and (21). Below we give, as examples, the explicit 



expressions for W 3 (t), Wi(t), and W 5 (t) 



W 3 (t) 
W 4 (t) 

W 5 (t) 



1 

3 
4! 
1 

4./,, 



5! 



dt'[W 1 (t'),[W 1 {t'),H 1 {t')]\, 



dt'[W 1 (t , ),[W 1 (t , ),[W 1 (t , ),H 1 (t') 



dt'[W 2 {t l ),{W l {t'),H l [t l ) 



df[Wi(f), 



[W X (*'), [W), (*')]]]] 

^[w^Miniwiinm')}}}. 



(27) 



As we show in the next section, the expansion obtained 
above may be useful in several cases and in particular for 



obtaining effective time independent hamiltonians, when 
the operator Unit') in the expansion can be approximated 
by the exponential of the product of the time with a 
constant operator. 

Notice that besides the fact that they are Hermitian, 
no restriction was made on the operators W n (t) for n = 
2, 3... until now. Special choices of W n (t), other than the 
one we have chosen to discuss in this paper, may lead to 
interesting applications in specific cases. 

III. EXAMPLES OF APPLICATIONS 

In this section we describe three examples of applica- 
tions of the method: i) the problem of a linear harmonic 
oscillator subjected to a driving force; ii) the Raman res- 
onant transition inside a cavity; iii) the ultrastrong cou- 
pling (USC) and deep strong coupling (DSC) regimes of 
the Jaynes-Cummings (JC) model. 

We start with the well known problem of a linear har- 
monic oscillator subject to a driving force —gf(t). The 
Hamiltonian is given by (h = 1) 



H = w(a f a + 1/2) +g/(t)(a + a f ) 



(28) 



where a and are the usual annhilation and creation 
operators satisfying the algebra [a,a^] = 1. 

We first take Hq = uj{a)a + 1/2) and go to the inter- 
action representation by defining U(t) = e~ lHot / h U~i(t), 
where 



dUijt) 
dt 



with 



Hi(t) = f{t){ae~ wt +tfe tult ). 



(29) 



(30) 



In this case, [H 1 (t), H^t')} = -2if(t)f(t') sinw(i - t 1 ) is 
a c-number. Therefore, W n = for n > 2, U 3 = I, and 

Wi(t)= f dt' f(t')(ae- iut ' +a}e iut '), (31) 
Jo 

W 2 (t) = i[ dh [ dt 2 /(ti)/(*a)sinw(ti-t2). 



o 







Then, the time evolution operator in the interaction pic- 
ture is given by 



U 1 (t) = e i ^D(v(t)) 



(32) 



where tp(t) = g 2 W 2 (t) is a time-dependent phase and 
D(v(t)) — e v ( t ) a1 ~ v (*) a is the displacement operator and 



v(t) = -ig / dt'f(t')e 
Jo 



(33) 



Another example is the case of resonant Raman scat- 
tering inside a cavity. Consider a three-level A atom in- 
teracting quasi-resonantly with a mode of frequency u>i 



4 



of the cavity field and a classical field of frequency u>2, as 
schematized in Fig (III). The two lower levels, \g) and |e), 



are closely spaced in energy and can make quasi-resonant 



dipole transitions to an upper level 



and 



uji e are 



the energy differences between the upper level and the 
lower levels |<?) and |e), respectively. The Hamiltonian 
that describes the interaction in the rotating-wave ap- 
proximation is given by 

H = H 'o+ H Lt 

H' = Ld ig \i)(i\ + Lo eg \e){e\ + uj{a)a 



H ln t — ^ifl 



{g\a - 



-iui 2 t 



i){e\ 



n ig \g)(i\tf +Sl ie e iu > t \e){i\ 



(34) 



where f2j 9 is the vacuum Rabi frequency associated to 
the cavity field of frequency ux, while O ie is the Rabi 
frequency associated to the external classical field of fre- 
quency LJ2- Assume that the initial cavity field state 
has a photon distribution with low photon average num- 
ber. In Ref . [5] , it has been shown that if the detuning 
5 = u>i g ~uji w uji g ~Lo eg —LU2 is such that \8\ ^> Q g i 3> Q, e i, 
it is then possible to show that the Raman transition 
\g, n + 1) f-> |e,no) is resonant for a certain n depend- 
ing on the detunings of the driving field. Here we rederive 
the conditions on the frequencies that make the process 
resonant and the effective hamiltonian for the system. 



i) 



9) 





e) 



FIG. 1. Raman transition of a A atom inside a cavity. 

Assume that the classical field frequency is tuned to 
W2 = Wl - UJeg - (n + l^/A + Sl 2 JA, (35) 
with A satisfying the equation 

A = cu ig -ux + + 2 ^)/ A + 2n O^/A, (36) 
where no is an integer and |A| 3> £l g i 3> Q e 



We first write the Hamiltonian of Eq. (34) as H 
Ho + Hi n t w i tn 



Ho — H o + H S s 

Hint = -ffi'nt — H S S, 



(37) 



where H' and H( nt are given by Eq. (34). Also, i?ss is 



given by 



ft 2 

- •• '";,/ 

A n 



Hss = { — nil 



A, 



^a\g)(g\ 



SI 2 
A, 



where 



g - wi - uj ig + fli g {2n + 1)/A 



A 

A e = LJ 2 — LU. 



nln/A. 



\e)(e\, 
(38) 



(39) 



We then write the time evolution operator in the interac- 
tion representation with respect to Hq and use our uni- 
tary perturbative expansion. Neglecting terms that vary 
very rapidly with time, we obtain 



XWx = -H ss t 
XW 2 = H S s+H eS t. 



where 



Therefore 



H e s = 



Uj(t) 



.M^Qe){ 9 \a+\g){e\tf) 



e +iH S st e -i(H ss + H cii )t 



(40) 



(41) 



(42) 



Using the Baker-Hausdorff formula and neglecting the 
term depending on the commutators of Hgg and H e g, we 
may write 



Uj(t,0) 



-iH c!! t 



(43) 



That is, -ffeff can be considered an effective Hamiltonian 
of the interaction picture associated to the choice of Hq 
given in the Eq. ( 37 1 . 
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FIG. 2. Survival probability of \g; 0} vs. r = cut. Solid line: 
exact solution; dashed line: U « UqUi; dot dashed line, first 
Born approx.; dotted line U « Uo- ijJq/co = 0.6;g/uj = 0.5. 

Consider now the situation of the Jaynes-Cummings 
model in the USC regime between a cavity mode and a 
qubit, g/uj > 0.1. This situation is currently accessible 
to experiments using superconducting qubits and cavities 
in circuit quantum electrodynamics |10Ullj . In this case, 
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FIG. 3. Probability of \g; 0} to make a transition to the state 
\e; 1) as a function of r = ut. Solid line: exact solution; dashed 
line: U as UoUi; dot dashed line: first Born approximation; 
dotted line U ~ Uq. uq/u = 0.6; g/u = 0.5. 



the rotating-wave approximation is no longer valid and 
one should consider the full interaction Hamiltonian 

H = uia/a + g(a j ' + a)a x + ujoct z /2 . (44) 

In the case where ujq = 0, it reduces to 

H' = wcfia + gia) + a)d x . (45) 

The eigenstates of H' are the product of displaced num- 
ber states [Hj and the eigenstates |±) of a x , associated 
with the eigenvalues ±1 



\±n;±)=D(Tx)\n)®\±), 



(46) 



where x — g/uj, D(v) — e va v a is the displacement 
operator, \n),n = 0,1,2.. are Fock states, and a x \±) = 
±|±). The eigenstates D(^fx)\n) of H' are degenerated 
and associated with the eigenvalue (nw — g 2 /u>). 

In basis {| ± n; ±), n = 0, 1, ...}, uj a z /2 is written as 

uj d z /2 = uj /2j2(n\D(2x)\m) \ + n; +)(-m; -| + H.c. 

n, m 

(47) 

In Ref. [13], it has been proposed an approximation 
which keeps only the terms with n = m'va the right hand 
side of Eq. (|47]), that is, 



Using the approximate Hamiltonian Hq as our zeroth- 
order approximation, we have found that our method is 
well suited for describing transition probabilities for a 
very large range of ojq/uj and g/u > 0.1, including both, 
the USC and DSC regimes of the JC model [H]. Let us 
write H = Hq + Hi nt , where 



H- mt = io a /2 (n\D(2x)\m) | + n; +){-m; -| + H.c. 

(50) 



From Eqs. (48) and (50) we can easily calculate 
\Hi(t) = e^H^te^^" 1 and the operators W x and W 2 



using the expressions given in Eqs. ( 10 ) and ( 23 ) 
1.0 f 
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FIG. 4. Survival probability of \g; 0} as a function of r = ojt. 
Solid line: exact solution; dashed line: U ^j- iS ot e -iXWiW . 

e - l H t e - l XW 1 {t) e -rXW 2 {t) dQtted rme 



dot dashed line: U 
U f 



lHot . ujq/uj = I; g/u = 0.; 
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H a = H' + uj a /2j2(n\D(2x)\n) \ + n; +)(-n; -| + H.c, 

n 

(48) 

where (n\D(2x)\n) — eT 2x ~ C n (Ax 2 ). The eigenstates of 
Hq can be easily written as 

-L(|+„; + )±|-n;->, (49) 

and are associated to the eigenvalues nui — g 2 /ui ± 
(LO /2)(n\D(2x)\n). 



FIG. 5. Transition probability from \g; 0} to |e; 1} vs. r = ut. 
Solid line: exact solution; dashed line: U w UoU\; dot dashed 
line: U ~ U0U1U2; dotted line U ~ Uq. uo/u = l;g/u = 0.8. 

In Figs. [2] and [3j we present curves corresponding to the 
exact result calculated numerically, the results calculated 
using three approximations for the time evolution opera- 
tor: e~ lHat , e - l #ot e -« Wi(t) ^ anc L Dyson approxima- 
tion to first order in ujq. The results show that both, the 
approximation H = H and the first Born approxima- 
tion, do not describe the transition probabilities for the 
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FIG. 6. Survival probability of \g;0) as a function of r = 
ut. Solid line: exact solution is superposed with U w UqIIi; 
dotted line U « Uq. wo/w = 0.5; g/o; = 2. 
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FIG. 7. Probability of \g;0) to make a transition to |e; 1} as 
a function of r = uit. Solid line: exact solution is superposed 
with U » UqUi; dotted line (7 » t/o- wo/w = 0.5; g/u> = 2. 

USC regime. On the other hand, our unitary expansion 
describes quite well the results, even when we take only 
the first non-trivial contribution. Note that the contri- 



bution of W2 is not shown and gives a small correction. 

In Figs. [4] and [5| we also present curves corre- 
sponding to the exact result calculated numerically and 
those calculated using three approximations for the 
time evolution operator: e~ iHot , e -»#o* e -iAWi(t) an( j 
e -i8 t e -i\Wi(t) e -i\w 2 (t)_ The firgt Born approx i ma tion 
is not shown and represents a small correction to the ap- 
proximation H s» Ho. We see that, although the terms 
associated with W\ gives the main contribution, the 
terms associated with W2 does improve the approxima- 
tion. These considerations are also applicable to the DSC 
regime [14] . as can be seen in Figs. 6 and 7. 

IV. CONCLUSIONS 

We introduced a perturbative unitary expansion for 
the time evolution operator as a product of exponentials 
of antihermitian operators. Consequently, this expansion 
can be truncated at any order of approximation while 
keeping unitarity. We have presented three examples: a 
harmonic oscillator with a time-dependent force, the Ra- 
man transition inside a resonant cavity, and the James- 
Cummings model in the USC and DSC regimes. 
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